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Abstract 

We consider a bistable integral equation which governs the stationary solutions of a 
convolution model of solid-solid phase transitions on a circle. We study the bifurcations 
of the set of the stationary solutions as the diffusion coefficient is varied to examine the 
transition from an infinite number of steady states to three for the continuum limit of 
the semi-discretised system. We show how the symmetry of the problem is responsible 
for the generation and stabilisation of equilibria and comment on the puzzling connection 
between continuity and stability that exists in this problem. 



1 Introduction 

Integro-differential equations are used to model various phenomena in materials science [U [21 
[21 [5l O [T7] and biology O [3 [211 [25| , which involve non-local diffusion/dispersal mechanisms. 
We consider the integro-differential equation (IDE) 

ut = £([ J^{x-y)u{y,t)dy-u{x,t) I J°^{x - y)dy) + f{u), (1) 

where the L^(]R) kernel J°° satisfies J°°(x) > 0, J°°{x) = J°°{—x) and /(n) is a bistable 
nonlinearity. Below we routinely consider f{u) = u{l — u^) and kernel 



J°°(x) = y ^ exp(-100x2), (2) 

so that J^J^ dx = 1. To obtain a well-defined problem, ([1]) has to be supplemented by 
a suitable initial condition, u(x,0) = uo{x) which needs to be chosen in a suitable function 
space, see [TS j [T9l [T6]. 

The convolution equation ^ is the L^-gradient flow of the free energy functional 

E{u) = \e f [ J°°ix-y) (n(y) - u{x)f dxdy + [ F{u)dx, (3) 
4 Jm. Jr Jr 

where F{u,t) is the smooth double well potential, F'{u) = —f{u). 
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For an overview of the use of ([T]) in materials science, see |10j . There are many papers deahng 
with the mathematical analysis of this equation, which examine existence and stability of trav- 
elling waves [3j, the structure of the stationary solutions set |2j, propagation of discontinuities 
pT] . coarsening [9j and long time behaviour [T5 l [T9l [Ml [20] , 

Note, in particular, that in [T3] it is shown that if the diffusion coefficient e is sufficiently 
large, a "Conway-Hopf-Smoller" type result holds: the only stable steady state solutions, 
say, in L°°(M), are the constant stable steady states of the kinetic equation ut = f{u). Thus, 
if we choose f{u) = u{l — u^), the stable states are u = 1 and u = —1. On the other hand, 
if e = 0, ([T]) admits an uncountable set of equilibria: let X, Y and Z be any disjoint sets 
such that A U B U C = M, then a function u{x) that is equal to 1 on X, —1 on Y and 
on Z is a steady state solution. Note that if Z = 0, all the resulting equilibria are stable 
in L°°{I). Furthermore, it is shown in [9] that there exists an > which depends on the 
kernel J°°, such that for all < e < Eq the set of steady state solutions of ([1]) is in one-to-one 
correspondence with the set of equilibria of ut = f{u). Hence, in view of the above, it is of 
interest to perform a bifurcation analysis of the set of steady states of ([TJ , 

= £ / J^{x - y){u{y) - u{x)) dy + f{u), (4) 



as we decrease e from some initially large value to zero, and investigate the transition from a 
finite to infinite set of solutions. 

To the best of our knowledge, such a study has not been performed before. The object of 
this paper is precisely such a study of the spatially discretised version of ([1]). For simplicity, 
here we restrict ourselves to 1-periodic patterns. 

If we choose spatially one-periodic initial data ti(x,0), then from ([T|) it is clear that for all 
a; € M and t G M+ 

u{x, t) = u{x + 1, t). 

Then from ([T]) we have 

ut = e J°°{x-y) {u{y, t) - u{x, t)) dy + f{u) 

oo „r+l 

= eY, J^{x-y){u{y,t)-u{x,t))dy + f{u) 

oo „i 

= £ ^ / J°°{x- z-r){u{z + r,t)-u{x,t))dz + f{u) 

r=-oo 

= £ / J{x — z) {u{z,t) — u{x,t)) dz + f{u), (5) 

^0 



where 



J(x) = J°°{x-r) (6) 



and x E [0, 1]. Thus, for 1-periodic initial data we only need to solve the problem ([T]) on the 
interval Vt = [0,1] with the kernel J{x). For the kernel given by ([2]), J°°{x) and J{x) are 
plotted in Figure [H 



2 



Lemma 1. For J defined by ^ the following two properties hold. 
1. If J°°{x) = J°°(-x) we have that 



J{x) = J{l-x). 



2. Qj{x)dx = jr^J^{x)dx. 



Property 1. above has an important influence on the spectrum of the matrix governing the 
semi-discretised version of ([T]) as we explain in the next section. From now on we work on 
[0, 1] and use the kernel J given in ([6]) 



(a) 



(b) 






Figure 1: The kernels J°°{x) in (a) and J{x) in (b) for the case of equation 



2 The semi-discretised system 

We discretise in space using piecewise-constant functions [9J and collocating at the uniformly 
spaced element mid-points, x = x-,i , j = 0, 1, 2, • • • , — 1. Setting u,- = uix-.i ,t), we have 
the semi-discrete approximation of ([T|) given by 



where now u{t) € M^^, supplemented with some initial condition m(0) = uq ^ 



(7) 



The 



nonlinearity F : R-^ i— > is given by Fj{u) = f{uj). It remains to specify the N x N matrix 
An- If we put h = j^, its elements are given by 



hJ{\xj^i\) 



h 



(8) 



From Lemma [T] it follows that jd^r is a symmetric circulant matrix generated by the elements 
ai^i, . . . , ai^N- Hence the theory of circulant matrices can be used to characterise its spectrum 
precisely. Let Wk be the distinct roots of - 1 = 0, so Wk = exp (^^), for k = 
0, 1, 2, N — 1. Then the following theorem holds: 
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Proposition 1 ([23]). Let be the circulant matrix defined by ai^i,ai^2T ■■■,0'i,N ■ Then 
—Ai\f is diagonalisable with eigenvalues 



Afc = -Ki + ai,2Wfc + ai,3VF| + ... + ai,NW^-\ (9) 

with corresponding eigenvectors = [l,Wk,W'^, ...^W^~^^ . 

Let us see what this imphes in our case for the spectrum of the discretisation. 

Lemma 2. The following three properties hold for the spectrum of A]\j 
1- Ao = 0; 

2. Afc = Xn-u; 

3. Let L]sf be the convex hull of the set of non-zero eigenvalues of A]\f. As N — )• oo, Ijv 
converges in the Hausdorff metric to the set 



Loo 



1 

J{x)dx 



I J{x) ex.p{2TTix) dx, / J{x)dx. 

JQ 



Before we prove this lemma, let us explain what it means. First of all, we must have a zero 
eigenvalue with a constant eigenvector, because, like in the case of the Neumann Laplacian, 
the equation 

poo 

ut= J°°{x - y){u{y,t) - u{x,t)) dx, 



conserves mass. 

Secondly, the pairing of the eigenvalues is simply the consequence of the symmetry J{x) = 
J(l — x) inherited from the evenness of the kernel J°°. Finally, the third part of the lemma 
implies that as N ^ oo, the spectrum accumulates at the point J{x) dx. Note that in the 
case of J°°{x) = y^IOO/vt exp(— lOOx^), we explicitly have 

loo = [1 - exp(-7rVlOO), 1] = [0.094, 1]. 

Proof. 1. From ([5]), putting Wq = 1, we immediately obtain from Q that Aq = 0. 

2. From part 1. of Lemma [1] it follows that for all j = 2, . . . , A^, 

so that the matrix Aj\f is symmetric. Hence its eigenvalues are real. But then taking 
complex conjugates of AjyVk = XkVki we get that A^vj: = XkVk, or in other words AnVN-k = 
\kVN-k and hence = XN-k- 

3. Finally, by taking the limit as — > cx) in ([9]) we immediately obtain that 



Afc 



/ J{x) dx — J{x) exp{2TTikx) dx, 
Jo Jo 



k = l,2,.... □ 

Our aim is to examine bifurcations in this system and, below, we perform a numerical path- 
following of solution branches. Some of these will, by symmetry, arise in pitchfork bifurcations 
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from the trivial solution u = 0. Here we examine analytically the values of e where such bi- 
furcations may occur in the semi-discrete system and later we can compare to the numerically 
found values. Linearising around the zero solution, we have the eigenvalue problem 

eAnv + gradF(0)t> = ^iv, (10) 

and hence bifurcations from the zero solution will only occur if /x = 0, or in other words, if 

—A]\fV = V. 

e 

Thus, for the semi-discrete system ([7]) we can fully characterize the values of e where bifur- 
cations of the zero solution occur, namely 

ek:=^, k = Q, ...,N -I. (11) 

For example, for = 32, J'^{x) = ^ ^ exp(— lOOa;^) and f{u) = u{l — u^), we have using 
(j9|) , the results of Lemma [2] and the formula (jlip that bifurcations from the zero solution are 
expected at the values of e as in Table [2j Note that for this case of = 32, the value of ei 
agrees to 12 decimal points with the limiting value of ei, 1/(1 — exp(— vr'^/lOO)) as N ^ oo, 
(see part 3 of Lemma [2]). 



£l = £31 


10.6403416996149 




£9 = £23 


1.00033746722800 


£2 = £30 


3.06584313146254 




£10 = £22 


1.00005172586163 


£3 = £29 


1.69885748860222 




£11 = £21 


1.00000650971543 


£4 = £28 


1.25968856776757 




£12 = £20 


1.00000067252291 


£5 = £27 


1.09266327932320 




£13 = £19 


1.00000005703325 


£6 = £26 


1.02948119722480 




£14 = £l8 


1.00000000397031 


£7 = £25 


1.00800141737791 




£15 = £17 


1.00000000022729 


£8 = £24 


1.00180943793728 




£l6 


1.00000000002128 



Table 1: For = 32 Bifurcation values in terms of e of the zero solution. 



Let us examine the eigenvectors of —Ai\f in some more detail. Since both and fAr_fc are 
eigenvectors, we immediately have that Re (f^) and Im {vk) are eigenvectors. Define the cyclic 
shift a on u = {ui, . . . , u^) S IR'^ by 

a{u) = {uN, Ui, UN-l), 

then we have 

Lemma 3. If v is a real eigenvector of —A^ corresponding to a double eigenvalue X, then so 
is cr{v). 

This follows since if v is an eigenvector, then so is e^'^'^/^v. 

Remark. In the above argument, we can pass to the limit as — )• 00 and arrive at the 
somewhat startling conclusion that cos(27rA;x) and all their translates are eigenfunctions of 
—A = — (/q^ J{x — y){u{y) — u{x))dy) no matter what the kernel J{x) is as long as it has 
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the right symmetry property. Of course, cosines are also the eigenfmictions of the Neumann 
Laplacian. It is very pleasing to obtain such a result via a semi-discretisation. 
Finally we note that fixed points of the semi-discrete problem satisfy 



= eANU + F{u). 



(12) 



Thus at e = stable solutions are given by 



u = 



{ 



1 X eX, 
1 xeY 



(13) 



where X U 1" = [0, 1]. Unstable solutions at e = are given by 



u = 



1 xeX, 

1 xeY, 

X £Z 



(14) 



where X UY Li Z = [0, 1] with some nonempty Z. 

We use this to define solutions with different numbers of interfaces. When X = [0, a) and 
Y = [a, 1], < a < 1, we call u a one-interface solution of (I12p if for e = for some 
n G [0, — 1], cr"(n) = ai on X = [0, a) and cr"(n) = 02 on y = [a, 1], oi 02 € { — 1,0, 1}, 
ai ^ a2- That is loosely speaking we have at e = one jump in the solution upto cyclic 
shift. Two-interface, three-interface solutions, etc., are defined similarly. Thus, for example, 
the branch of solutions corresponding to orbit A in Table [2] are of one-interface and those 
corresponding to E are of three-interface. 

3 Results 

We take for our computations the kernel function 



with /(u) = u{l — u'^) and vary the parameter e. For small values of N it is possible to 
enumerate all possible solutions of the semi-discrete system ()12p with e = and to analyse 
their continuation to e > using the theory of bifurcation with symmetry. This we do below 
for = 4 and these analytic results were used to check the validity of our numerics. 
We implemented in Matlab a standard pseudo arc-length continuation algorithm with step 
size control as described in [I4l[22l[l3] for the discrete problem (fT2]) . Since An is a circulant 
matrix, we take advantage of reducing storage costs as the full information oi Aj\f can be 
obtained storing one row or column only, see |4] and references therein. Furthermore the use 
of the FFT for each matrix vector multiplication reduces the computational cost. We detect 
bifurcation points by observing where eigenvalues of the Jacobian J = D^F of the nonlinear 
system F(u, /3) = cross the imaginary axis and perform branch switching at those points 
by perturbing in the direction of the associated eigenvector. 
The arc-length £ of u{x) G C^(R) is defined in the standard way 
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and we approximate the arc-length of u{x) with the mid-point rule and using the standard 
forward difference approximation for the derivative. With a uniform discretization we get 




(15) 



j=0 



where h = Xj+i — xj and Uj ~ (xj). Note that although £ only makes sense for u £ 
however we can evaluate ih even when u is discontinuous at grid points. 
Then, for = 32 we compute the bifurcation diagram numerically and gain insight into the 
structure of the bifurcation diagram of the original continuous problem. 
Finally, we examine the large N limit and formulate the results of the numerics as two 
conjectures concerning the interplay of continuity and stability and the behaviour of saddle- 
node bifurcations as a —t- 1/2. 

For the continuous system, the symmetry group is 0(2) x Z2, and so for a finite number of 
nodes A^, we use Fjv = D^^ x Z2 equivariance structure [121 ITS]. 

3.1 The = 4 case 

If = 4, there are a total of 81 possible steady states at e = 0, 16 of them stable. The 
group r4 is generated by the shift p, the flip / and the reversal m. In other words, if 
V = (fi, V2, V3, V4), we have that 



Inverses of the nonzero eigenvalues of the 4x4 matrix A4 are {91.82, 183.63, 183.63}, so we 
expect primary branches to bifurcate from the zero solution at those values of e. Note that 
all primary branches have zero mean, but the converse is not true. 

Since here we know all the solutions at e = and their stability, and since symmetry properties 
are conserved on primary branches, we can cut down the work considerably by looking only 
at orbits of solutions under T4. In the table [21 we collect all the orbits, their lengths and 
the corresponding isotropy subgroups Ti^. There, (/) stands for the group generated by 



Now we can immediately draw the bifurcation diagram using the following three rules [121 11] • 
First a bifurcating branch must have the isotropy subgroup which is a subgroup of the isotropy 
subgroup of the primary branch; secondly dimensions of unstable manifolds have to match at 
a bifurcation point to satisfy the principle of exchange of stability, and thirdly at e = 0; the 
number of nodal domains must increase from one bifurcation point to the next. 
With these rules there is only one way to construct the bifurcation diagram; see Figure [2] 
(a) and (b), where the y-axis is not to any scale, and is only intended to make clear the 
end-points of various branches at e = 0. These figures show the bifurcation structure arising 
from bifurcations of the zero solution. 

We would like to make the following observations. The stable non-zero-mean branches cor- 
responding to the orbit K have to arise through a saddle-node bifurcation. Numerically, 
this happens at a value of epsilon ~ 49.294 that is smaller than the value £4 = 122.432 at 
which the branches of the orbit A become stable, see Figure [3] which shows the numerically 



p{v) 

fiv) 
m{v) 



{V2, f3, Vi, vi); 

: {-Vi, -V2, -Vs, -W4); 



/ G L>4 X Z2. 
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Name 


Orbit 


length 






(0,0,0,0) 


1 


L»4 X Z2 




(1,1,1,1) 


2 




A 


(-1,1,1,-1) 


4 


{f,p^m) 


B 


(1,0,-1,0) 


4 


{mp^Jp) 


C 


(0,0,1,1) 


8 


if) 


D 


(0,-1,1,1) 


16 


(I) 



Name 


Orbit 


length 




E 


(-1,1,-1,1) 


2 


{pm, fp) 


F 


(0,1,0,0) 


8 


ifp) 


G 


(1,0,0,-1) 


8 


(mf) 


H 


(0,1,0,1) 


4 




I 


(0,1,-1,1) 


8 


ifp) 


J 


(0,1,1,1) 


8 


ifp) 


K 


(-1,1,1,1) 


8 


ifp) 



Table 2: Steady states for = 4, the length of the orbits and isotropy subgroups T^x- We 
have separated the solutions into the homogeneous states, those connected with the first and 
second bifurcation of tt = and the two solutions connected by a saddle-node J and K. See 
also Figure [21 



(a) 



'x,^ 1 
3 



■.\ 1 



183.63 



E 
I 

H 
G 

F 



(b) 




2 



3 " 



\3 

4 \ 3 



91.82 e 



Figure 2: Bifurcations from the zero solution from the first (a) and second (b) bifurcation 
points. Solid lines represent stable and broken lines unstable unstable solutions respectively. 
We indicate above each branch the dimension of the unstable manifold. We do not show 
the homogeneous solutions or the orbits J and K (which are connected by a saddle-node 
bifurcation) . 
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Figure 3: Numerically computed bifurcation diagram from the zero solution from the first 
(a) and second (b) bifurcation points. Solid lines represent stable and broken lines unstable 
unstable solutions respectively. We indicate above each branch the dimension of the unstable 
manifold. We show the homogeneous solutions or the orbits J and K connected by a saddle- 
node bifurcation on (b). Compare to the theoretical prediction in Figure [2j 

computed diagram. We will see the equivalents of these statements in higher dimensional 
discretisations. 

Finally, we did not perform a Liapunov-Schmidt calculation to determine the order of bifur- 
cations at the double eigenvalue point e = 183.63, but the opposite assignment of stabilities 
cannot be reconciled with the above rules of bifurcation. 

3.2 The case of = 32 

Though an analysis similar to that in the case of = 4 can be attempted here, the numbers 
of orbits are astronomical, and we rely on our numerical continuation method, the results of 
which match exactly the predictions of the analysis in the case N = A. In Figure |4] we plot 
in (a)-(d) sample solution branches of the bifurcation diagram with N = 32. If we start with 
a large value of e we see in (a) and (b) the first bifurcation arises at e w 10.64 as predicted 
by the theory in Table [2j In (a) we show the continuation of the zero mean one-interface 
which undergoes a pitchfork bifurcation at e ~ 2.012. In (b) we have plotted the one-, three-, 
five- and seven interfaces and their stabilization. In (c) we show details of the bifurcation 
structure close to the pitchfork at e ~ 2.012 (note for clarity one branch of the pitchfork 
seen in (a) is not plotted). As a — )■ 0.5, the saddle-node bifurcation points converge to £3^2- 
This structure is repeated for the other n-interface solutions and is illustrated in (d) for the 
three-interfaces solutions. Here we see that the zero-mean one-interface solution branches 
stabilize at £3^2 ~ 2.012. Below we will formulate a conjecture concerning the limiting value 
which we call e? at which the one-interface branch with zero-mean stabilizes as A^ ^ 00. 
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4 The limiting problem and conclusions 



It is not hard to prove (see for example [3J) that if e > 1, steady state solutions of (jH) are 
continuous, since the function —£u + f{u) is monotone. Hence it is interesting to understand 
when the solutions lose continuity (certainly, for e = there are no non-constant continuous 
solutions) . 

The non-trivial stable one-interface zero-mean solution branches {a = 0.5) that originate at 
£ = were investigated in detail as we change N. If we define M by 



M := max 

j 



h 



then for a function this converges to max^.g[o,i] Wx\ and so we can identify where the 
solution is continuous. 

Figure [5] plots in (a) M against £ along a branch of one-interface zero-mean solutions for 
N = 2P , p = 4, 5, 6, 7, 8, 9, 10. If we let be the value of e at which this branch of solutions 
becomes discontinuous, then this figure suggests ef = 1. This is supported in (b) which shows 
for different £ convergence of the derivative M with on a log log scale. 
Furthermore the loss of continuity appears to coincide with the loss of stability. In Figure [6] 
we show numerically that the bifurcation values converge to = = e|] = 1 as iV — )• oo. 
Now we can collect our observations and form two conjectures. First we consider the zero- 
mean interface branches. Let be the value at which the zero-mean 2m — 1-interface solution 
becomes stable. Let be the value of £ at which this branch becomes discontinuous. Then 
we have 

Conjecture 1: = ej^. 

We can prove a very weak form of this conjecture for m = 1. From the results of [3j it follows 
that discontinuous stationary solutions will exist for any e such that the function 



g{u) := -£u 



[' J{s)ds + f{u) 
Jo 



is non-monotone. On the other hand, from Theorem 2.1 of [6j it follows if g{u) is monotone, 
there are no nonconstant minimizers of the energy functional ([3]). Hence we have ef < ef. 
However we do not have the inequality the other way. 

We now consider the saddle-node bifurcation of the non-zero mean interface solutions. Now, 
let Us be a branch of 2m — 1-interface stable solutions of ^ with mean s, and let £^ be the 
value of e at which the saddle-node bifurcation giving rise to the branch occurs. Then we 
have 

Conjecture 2: lim e^* = e^L. 

s->0 

These two conjectures, if true, would lead to the bifurcation picture sketched in Figure [71 In 
(a) we plot the zero-mean one-interface branch and have indicated the continuum of saddle- 
node bifurcations ef that approach the bifurcation at ef = ef = 1. In (b) we indicate the 
first four branches of the infinite number that bifurcate from zero, the branches of associated 
saddle-node bifurcations and here we have that liuis^o^ni = — ^m- addition our 
numerical investigation seems to indicate that ef = ef = 1 = e^ = e2. 

Finally let us consider the stable solutions - that is the solutions we expect to see from any 
simulation. Thus we have for e > 1 two stable solutions, then a region of parameter space 
with an infinite number of stable solutions of one and three interface type, then a region of 



10 



parameter space with one, three and five interfaces and so on. In conclusion the diffusion 
coefficient s determines the number and type of stable solutions. 
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Figure 4: (a) Continuation of the zero-mean one-interface (a = 0.5) solution branch. When 
£ ~ 2.012 shown by 'o' there is a stabilizing pitchfork bifurcation, (b) Stabilization of one-, 
three-, five-, and seven-interface solutions, (c) This is a blow up of the bifurcation diagram 
around the stabilizing bifurcation at e ~ 2.012 shown by 'o' for the a = 0.5 curve. Here stable 
solutions are continued from e = with different ratios (a) of —1 and 1 values of a. (d) A 
similar structure is observed starting from e = with two-interface solutions with different 
ratios (a) of —1 and 1. 
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(a) (b) 




Figure 5: In (a) we show the maximum of \ux\ in [0, 1] of the solution u{x) for one-interface 
initial data with ±1 and a = 0.5. In (b) we see maxx{ux) for different system sizes for some 
particular e's with system size as the x-axis with a clear change in behaviour at e = 1. 




Figure 6: (a) Pitchfork bifurcations of the first (2m — l)-interface solutions as N increases, 
(b) A log-log plot of e — z versus showing convergence to z = 1, 1, 0.9, 0.7839 for the one, 
three, five and seven-interface solutions. 
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Figure 7: (a) Proposed bifurcation diagram for the one-interface solutions in the N = oo 

case. We show the primary branch and the continuum of saddle-node from the saddle-node 
bifurcations. In (b) we indicate that this structure is then repeated for the three-interface, 
five- interface,... solutions. The solid circles represent lims^oSm' = — ^m- The bifurcations 
for the one and three interface solutions both occur at e = 1. 



15 



